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Abstract: 

Borexino, a liquid scintillator detector at LNGS, is designed for the detection of neutrinos and an- 
tineutrinos from the Sun, supernovae, nuclear reactors, and the Earth. The feeble nature of these 
signals requires a strong suppression of backgrounds below a few MeV. Very low intrinsic radio- 
genic contamination of all detector components needs to be accompanied by the efficient identifi- 
cation of muons and of muon-induced backgrounds. Muons produce unstable nuclei by spallation 
processes along their trajectory through the detector whose decays can mimic the expected signals; 
for isotopes with half-lives longer than a few seconds, the dead time induced by a muon-related 
veto becomes unacceptably long, unless its application can be restricted to a sub-volume along 
the muon track. Consequently, not only the identification of muons with very high efficiency but 
also a precise reconstruction of their tracks is of primary importance for the physics program of 
the experiment. The Borexino inner detector is surrounded by an outer water-Cerenkov detector 
that plays a fundamental role in accomplishing this task. The detector design principles and their 
implementation are described. The strategies adopted to identify muons are reviewed and their 
efficiency is evaluated. The overall muon veto efficiency is found to be 99.992 % or better. Ad- 
hoc track reconstruction algorithms developed are presented. Their performance is tested against 
muon events of known direction such as those from the CNGS neutrino beam, test tracks available 
from a dedicated External Muon Tracker and cosmic muons whose angular distribution reflects the 
local overburden profile. The achieved angular resolution is ^3° -5° and the lateral resolution is 
^35-50 cm, depending on the impact parameter of the crossing muon. The methods implemented 
to efficiently tag cosmogenic neutrons are also presented. 

Keywords: Cerenkov detectors, Neutron detectors, Particle identification methods, Particle 
tracking detectors, Particle detectors, Front-end electronics for detector readout, Large detector 
systems for particle and astroparticle physics. 
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1. Introduction 

Borexino is a large-volume liquid- scintillator detector for low-energy neutrinos and anti-neutrinos 
coming from the Sun, from supernovas, from nuclear reactors, and from the Earth (geo-neutrinos). 
The very low radioactive background conditions achieved for the scintillator allowed us to widen 
the scope of the experiment beyond the measurement of the monoenergetic solar 7 Be neutrinos flU]: 
measurements of solar 8 B neutrinos [Q] and of geoneutrinos [0] have recently been published. It is 
planned in the future to look also for solar pp, pep, and CNO neutrinos. 

A key element to reach the required sensitivity for these low-rate and low-energy event searches 
is the suppression of the background events created by natural radioactivity and cosmic muon ra- 
diation, so muon identification and tracking are therefore crucial. 

The Borexino detector was designed to have very low intrinsic background. The central scin- 
tillation volume, 278 1 of ultra-pure PC (pseudocumene) doped with 1.5 g/1 of the fluor PPO (2,5- 
diphenyloxazole), is contained in a spherical Inner Vessel (IV), 8.5 m in diameter, made of 125 jam 
thick nylon. It is floating in two buffer layers consisting of PC and a small amount of the light 
quencher DMP (dimethylphthalate). These inactive shells provide shielding from long-range y- 
rays and Rn emanating from external parts of the detector. The surrounding Stainless Steel Sphere 
(SSS) of 13.7 m diameter holds 2212 inward-facing 8" photomultiplier tubes (PMT) that detect the 
scintillation light from the central region. All these components form the Inner Detector (ID) [gj]. 
Due to the radiopurity of all materials and the self- shielding capability of the liquid scintillator, 
background levels corresponding to rates of less than one count per ton-day have been achieved in 
the central part of the IV (~ 100 1), the Fiducial Volume (FV). 

A second central requirement is the rejection of cosmic muons crossing the detector and the 
subsequent signals caused by muon spallation products. Though Borexino is located deeply under- 
ground, in Hall C of the Gran Sasso Laboratory (LNGS) 1 , the residual muon flux is ~ 1.2|i/m 2 /h, 
still too large for neutrino measurement, so the muons must be individually tagged. To accomplish 
this task the ID is surrounded by a powerful muon detector. This is composed by a high domed steel 
tank of 18 m diameter and 16.9 m height filled with 2 100 1 of ultra-pure water and instrumented 
with 208 PMTs which detect the muon Cerenkov emission. The Water Tank (WT) also serves as 
additional passive shielding against external radiation. This system is called the Outer Detector 
(OD). Moreover, the characteristic light patterns created by muons in both water and scintillator 
can be exploited to reconstruct the muon tracks. This is mandatory for the rejection of long-lived 
cosmogenic radioisotopes, such as n C, that constitute the main background for the forthcoming 
pep and CNO neutrino detection. Moreover, it is an important tool for the study of cosmogenic 
backgrounds. 

lr The rock mass above the detector corresponds to about 3 800 meters of water equivalent (mwe) and provides a 
reduction of the surface muon flux by about 6 orders of magnitude. 
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In this paper we describe the OD hardware and the identification and reconstruction of muons 
and of spallation neutrons. Sec. @ gives an overview of the components of the detector and of the 
electronics, the data acquisition system being described in [0], The criteria applied for the rejection 
of muons in both OD and ID are described in sec. |3|. The efficiencies of the individual methods, as 
well as the overall veto efficiency, are derived from the study of cosmic muons and of muons cre- 
ated by the CNGS (CERN Neutrinos to Gran Sasso) neutrino beam SJ. In sec. f| we present three 
track reconstruction algorithms we have developed for crossing muons. We evaluate their perfor- 
mance using muon tracks of known direction such as CNGS events observed also by the OPERA 
experiment, tracks registered by a small external auxiliary system we introduced for this purpose 
and cosmic muon and neutron samples. Finally, the detection of events in the wake of the passing 
muons is discussed in sec. |5] which is dedicated to cosmogenic neutrons. As the light generated 
by muon-induced particle showers severely interferes with the ID electronics, special techniques 
have been developed to efficiently detect both neutrons and subsequent decays of cosmogenic ra- 
dioisotopes. A brief conclusion is given in sec. |5|. The CNGS neutrino beam and the selection of 
corresponding events in Borexino are described in app. 

2. The Outer Detector 

The OD forms the outermost layer of the Borexino detector. PMTs distributed on both the SSS 
outer surface and the OD floor register the light produced by passing muons (see fig. |T|). 

This section describes the hardware of the OD, in particular the reflective foils enhancing the 
light collection, the physical arrangement and encapsulation of the PMTs, their LED calibration 
system, the electronics read-out chain, and the digital trigger system. 

2.1 Tyvek reflector 

The steel surfaces of both the SSS and WT are covered with reflective Tyvek 2 sheets in order to 
enhance the light collection (fig. @). The foils are a white paper-like material, ^200 jam in thick- 
ness, composed of pressed PET fibers (PolyEthylene Terephthalate). The spectral region relevant 
for Cerenkov light detection is defined both by the initial Cerenkov spectrum and the quantum ef- 
ficiency of the PMTs (fig. 0). In this region the reflectivity of Tyvek is twice that of steel (fig. |]). 
Convoluting the spectrum of reflectivities with the PMT quantum efficiency and integrating returns 
an effective reflectivity of 38 % in the case of steel and 79.5 % for the Tyvek foils. Although there 
is a clear gain in total light collection, and therefore in veto efficiency, these foils make the muon 
track reconstruction in the OD more complex because of stray light (sec. |4.1.1|) . 

2.2 Photomultiplier tubes 

The OD is equipped with 208 8" ETL-9351 PMTs manufactured by Electron Tubes Limited (for- 
mer Thorn Emi ETL) fl7|]. Nominal operation values are listed in tab. |T|. The wavelength-dependent 
quantum efficiency is presented in fig. where the Cerenkov light spectrum is also given for com- 
parison. 

2 Registered trademark of the DuPont company. 
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Figure 1. The Borexino Outer Detector. PMT positions are shown, grouped in rings. In the table, ring 
positions are indicated as stainless steel sphere (SSS) — Upper (UH) and Lower (LH) Hemisphere, the OD 
floor and the 45° -slope bordering the floor. The vertical z coordinate corresponds to the Borexino coordinate 
system in which the origin is in the center of the SSS. The last column gives the number of PMTs per ring. 
PMT belonging to the same ring are connected by black lines in the figure (sec. 





Figure 2. OD PMTs. The encapsulated photocathode is emerging from the Tyvek sheets. The calibration 
optical fiber is also visible. 



As the PMTs are operated in the aqueous environment of the OD, they are fully encapsulated 
to shield their base electronics from water. Fig. |3] shows the PMT contained in a conically shaped 




wavelength [nm] wavelength [nm] 



Figure 3. The 8" ETL-9351 PMT quantum effi- 
ciency and the Cerenkov light emission spectrum. 




Figure 5. Encapsulation and support structure of the 
OD PMTs. 



Figure 4. Tyvek reflectivity as a function of wave- 
length compared to steel. 



Table 1. Manufacturer specification for 8" ETL- 
9351 PMTs, operated at an amplification gain of 10 7 . 
SPE stays for single photoelectron. 



voltage 


-1650 V 


SPE rise time 


4 ns 


SPE FWHM 


7 ns 


SPE fall time 


8 ns 


linearity peak current 


10 mA 


SPE peak-to- valley ratio 


2.5 


photocathode sensitivity (420 nm) 


26.5 % 


transit time spread (FWHM) 


2.8 ns 


pre-pulsing (2a) 


3% 


after-pulsing (0.05-12.4 jus) 


2-5% 


dark counts (room temp.) 


3 kHz 



steel case: a transparent (^90 %) PET foil at the top allows the light to reach the photocathode, 
retaining an acceptance angle close to 180°. The foil is flexible so that the water pressure is passed 
on to the PMT. The PMT head and the dynodes structure are shielded by a conical |i-metal foil 
from Earth's magnetic field. The remaining space between the PMT and its steel case is filled 
with a transparent mineral oil, whose refraction index of n « 1.5 matches both the PET foil and 
the PMT glass. The only (unavoidable) discontinuity in the light transport arises between water 
(n = 1.33) and PET foil. The PMT socket is plugged into the voltage divider that is in turn cast 
in polyurethane and in this way firmly connected to the housing bottom. An additional clamp (not 
shown in fig. ||) holds the PMT in its place 3 . Both HV and signal are fed via a Jupiter connector and 
a 55 m long submarine HV cable to the outside counting room. PMTs mounted in the WT can be 

3 This is needed, in case the PMT is positioned upward-looking, in order to counter the buoyancy force that otherwise 
pulls it out of the socket. 
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seen in fig. g. The design and realization of the encapsulation proved to be effective, a fact testified 
by the contained failure rate: after loosing ^2.5 % of PMTs during the initial filling of the WT, for 
the first three years of operation the failure rate has been only ^2.5 %/year. 

While the original detector design foresaw an equal distribution of the 208 PMTs on the SSS 
outer surface, rearranging the lower fourth of the PMTs to the floor of the WT significantly en- 
hances the direct light collection for vertical muons and consequently the tracking possibilities of 
the detector. Monte Carlo (MC) simulations indicated that the veto efficiency was unaffected by 
the rearrangement. The final PMT deployment is described in fig. |T|: 154 radially outward-looking 
PMTs on the SSS are arrayed in eight horizontal rings on the upper 2/3 of the sphere; 34 PMTs are 
ordered in four concentric rings on the WT floor, while 20 PMTs are mounted on a steel slope of 
45° inward inclination that runs along the bottom of the domed steel tank; the slope (about 1 m in 
height) was introduced for engineering reasons. 

2.3 Electronics 

Borexino PMTs produce a single photoelectron (SPE) signal about 12 mV high and 7 ns wide 
(FWHM), corresponding to ^1.6 pC charge. However, in the OD, a single PMT can easily see 
up to ^100 PE if directly hit by the Cerenkov cone. An ad-hoc data acquisition chain has been 
designed to ensure the necessary dynamic range, without sacrificing too much the SPE resolu- 
tion. The system is based on Charge-to-Time Converters (QTCs) that integrate the PMT signals. 
The logic pulses whose duration is proportional to the input signals' area are read by commercial 
multi-hit Time-to-Digital Converters (TDCs). DAQ occurs through the VME dataway bus and is 
mastered by a Motorola PowerPC (PPC) computer running custom software [0]. 

2.3.1 High Voltage supply and signal decoupling 

All PMTs are operated at a nominal gain of 10 7 . The negative operating voltage required to achieve 
this gain was measured for each PMT in a dark room test facility before the installation in the 
WT. Typical values lie between ~ 1150V and ^1750 V and the current drawn at nominal gain 
averages ^70 |uA. To supply this power we use a CAEN 4 S Y527 mainframe with 9 CAEN A932AP 
distributor boards, with 24 channels each. The PMTs are AC coupled, i.e. the signal and the 
HV travel along the same cable. The high voltage decoupling is performed by custom passive 
decouplers that also house the distribution system of calibration pulses to all channels. After the 
decoupling stage, the signal features a bipolar shape, however the overshoot is very small and long 
(^ms) and its interference in the signal integration can be easily kept to a minimum. 

2.3.2 Front-end electronics 

The front-end electronics is made of 14 QTCs, charge-to-time converter custom boards. The board 
(fig. ^) is a 9U single VME unit which takes only power from the VME backplane with no connec- 
tion to the dataway bus. For each of its 16 channels the board issues a Primary Output (PO), a logic 
signal whose leading edge provides timing information and whose length (distance to trailing edge) 
is proportional to the charge. The board also has a Secondary Output (SO), a logic step function 
pulse whose height is proportional to the number of channels firing in coincidence. SO signals are 

4 CAEN S.p.A., Viareggio (LU), Italy; www.caen.it 



-6- 



Figure 6. One of the Outer Detector front-end QTC boards. 



sent to the Muon Trigger Board (MTB) for the Outer Muon Trigger (OMT) formation (sec. |2.3.4Q . 
Both the width and the height of the steps can be trimmer-adjusted on the boards. Finally, the 
board has a Tertiary Output (TO) issuing an analog sum of the input signals with a built-in fixed 
amplification factor of about two 5 . 

The input signal is split in three paths reaching the QTC inputs (amplif.: x5), the threshold 
comparator (amplif.: x 1 1) and the TO summing block. 

The heart of the channel block is the LeCroy 6 MQT300 QTC chip. The comparator is operated 
with a threshold of -20 mV, corresponding to ^0.2 PE; its logic pulse gates the QTC for a time At, 
set to ~ 110 ns, a value optimized to preserve linearity over the whole dynamic range. The chip 
integrates the input pulse together with Vo, its offset (or pedestal), issuing a logic pulse whose 
length is given by: 



where P is the pedestal length in ns and G is the channel gain, trimmer-adjustable in the range 
[10-60] ns/pC. We use the lowest possible gain that ensures a dynamic range beyond 100 PE (cor- 
responding to ~ 2us pulse). Similarly a trimmer allows tuning of the pedestal length: typical values 
are ^700 ns. The comparator is disabled during the PO issuing and re-enabled immediately after. 
The dead time of the channel is therefore the time occurring between the closure of the integration 
gate and the end of the PO signal. 

The QTC boards are daisy-chained by the control bus running along a single ribbon cable. 
The MTB is the bus master and its behavior is fully programmable through the PPC. The control 
bus includes an ENABLE line, a test pulse, controls setting the run mode and integrity signals 
(described in sec. |2.3.3|) . 

5 These outputs were introduced as an additional analog trigger system, however the performance of the digital trigger 
proved to be beyond expectations making the exploitation of the TO unnecessary. 
6 LeCroy Corporation, Chestnut Ridge, NY (USA); www.lecroy.com 
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Two special run modes, enabled via control lines, allow the gate width and pedestal length to 
be sent through the PO. This is done at the beginning of every run. The pedestals are then subtracted 
during data processing and the gate widths are available for precise charge calibration. 

2.3.3 Digital electronics 

The digital electronics of the OD is made of two CAEN v673 TDC units. The unit houses 4 TDC 
chips, each one with 32 channels, yielding an overall 128 channels, with a time resolution of LSB 
~ 1.04 ns (Least Significant Bit). 

The TDCs are multi-hit devices with the ability to record up to 16 signal edges (rising or 
falling) for each event with a time depth programmable up to 64|as. We select 16 |us to closely 
match the trigger gate to the ID electronics 0]. The minimum time distance between 2 successive 
edges is 8 ns, a requirement well met by our QTC signals, typically in the range ^[0.7,2] jus. 

The unit keeps digitizing incoming pulses but data are periodically overwritten unless a trigger 
occurs. The unit works in COMMON STOP mode, with the trigger signal (sec. |2.3.4|) , properly 
delayed, being the STOP. Each chip has 4 event memory buffers that allow deadtime-less operation 
but which are easily filled by a fast sequence of triggers. STOP signals arriving when the 4 buffers 
are full are ignored by the board. For this reasons the ability of the DAQ to empty the buffers in 
the shortest possible time is particularly critical. 

The TDC control bus is daisy-chained among the boards, mastered by the MTB (sec. |2.3.4Q 
and includes STOP and RESET signals to TDCs as well as FULL and BUSY flags from the TDCs. 
When the BUSY condition is on either board, the MTB issues a VETO to the main trigger to avoid 
DAQ jamming. 

Synchronism among the 8 chips is not guaranteed by the units themselves: a FULL or a BUSY 
condition may regard only one or a few chips, while the others keep recording data. Moreover a 
read-out operation may result in data chunks belonging to different events according to the presence 
or absence of data in a chip for a given event. In our design the synchronism is insured externally 
by the presence of integrity channels. These are the first channel of each QTC board, pulsed in 
time with the trigger signal to guarantee a minimum of one hit per event in every chip. In this way 
event buffer occupancy of the 8 chips can be described by a single variable. 

Read and write operations are independent. The first word is the value of a counter holding the 
number of events occurred since the last read operation. Following hit words contain the difference 
of the gray counter readings upon each hit and common stop arrival, i.e. the negative hit time 
relative to the trigger signal in clock ticks. The board supports sequential as well as block read-out 
as described in [0], 

2.3.4 Trigger system 

In Borexino the trigger condition is a general condition for the whole system (Borexino General 
Trigger, BGT) which means that the whole detector is acquired simultaneously, no matter whether 
the trigger condition was risen by inner or outer detector. 

The OD features a trigger subsystem, fully integrated in the main one which receives and 
handles the BGT and that generates the OMT signal and forwards it to the main trigger system. 
These tasks are carried out by the MTB, a custom 6U VME unit (fig. 0). The board handles both 
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Figure 7. The Muon Trigger Board (MTB), before front panel application. 

QTC and TDC control busses (sec. |2.3.2| and sec. |2.3.3|) . The lines can be set or checked through 
VME by the PPC software. 

The board interface receives from the main trigger system the BGT and the TRGID, a 16-bit 
pattern to flag data with a universal event number. The board in turn sends out the OMT and the 
OR of the two TDC BUSY flags. 

Upon BGT arrival, the MTB issues the integrity signal to the QTCs, the COMMON STOP to 
the TDCs (delayed ^70 ns with respect to the integrity signal), acquires the TRGID pattern and 
finally raises a VME interrupt, within 200 ns of the COMMON STOP pulse, informing the DAQ 
that an event is ready to be read. 

The task of OMT formation is delegated to a daughter board for noise isolation. As described 
in sec. |2.3.2| , the QTC boards issue (SO) a step function encoding the number of channels fired in 
the board. The width of the steps defining the trigger time window is usually set to 150 ns, while the 
height is set to lOOmV/channel. Signals are summed and discriminated against a programmable 
threshold (8-bit DAC), which can be set at run start through the run controller. The inputs have a 
dynamic range that spans up to +2 V before saturation occurs, allowing the threshold to be set up 
to almost 20 channels. During normal data taking we use the equivalent of 6-7 PMTs (see also 
sec.gXTl). 

2.4 Calibration System 

The calibration system of the PMTs must provide both a relative time alignment with sub-ns preci- 
sion and a charge calibration on SPE level. The accuracy in time calibration has a large impact on 
the performance of the OD muon tracking (sec. |4.1.1|) , but has also influence on the efficiency of 
the muon identification flags (sec. PTT1) . Providing a synchronous time-signal on all photocathodes, 
relative offsets and jitter of both the PMT and the subsequent electronics chain can be determined. 
The calibration of the PMT SPE response ensures the comparability of photon numbers detected 
in individual PMTs and is necessary for the determination of the PMT center-of-charge that is the 
basis of the OD tracking. 

The OD calibration system is based on external LEDs. The light of individual LEDs is fed 
by 55 m long optical fibers 7 to the front of the OD PMTs (fig. |). The UV-LEDs 8 feature a peak 

7 Huber-Suhner (Germany), commercial fiber, length: 55 m, diameter: 200 (am, SMA connectors. 
8 Roithner Lasertechnik (Vienna, Austria), NSHU550 LED 
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emission at 370 nm and a time jitter of 0.8 ns. They are mounted on 18 customized VME boards 
of 12 LEDs each. The boards include programmable delay lines on every channel to correct inter- 
channel time offsets and allow for the possibility to adjust individually the voltage across LEDs. 
The individual LEDs are cross-calibrated using an outside PMT of the same type as the OD PMTs. 
This makes it possible to identify time and voltage parameters that achieve a synchronized light 
signal of equal intensity. This operation occurred immediately before the beginning of the data 
taking and could be repeated upon necessity during a maintenance shutdown. 

A custom control software accessed via the run controller allows the timing and intensities to 
be varied during DAQ if necessary. A weekly calibration of the OD PMTs is performed, monitoring 
drifts in the timing and charge of SPE signals. 

Fig. |8]a shows the sum signal of 10 5 LED calibration events, corresponding to a typical cal- 
ibration run. The solid red (black dashed) line shows the overall time distribution after (before) 
applying the time correction computed from the individual channels offset from the overall mean. 
Fitting the main peak with a Gaussian, the width improves from 3.5 to 2.1 ns (lcr). This directly 
corresponds to the improvement in relative PMT timing. 

Fig. |8|b shows the SPE charge spectra before (black dashed) and after (solid red) calibration. 
The first was corrected by the mean charge of all channels to allow a comparison of the distribu- 
tions. Here, the calibration gives correction factors of 0.5 to 2 to equalize the average SPE charge 
to a common value. As a consequence, the SPE peak is more prominent after calibration. 




760 780 800 

time to reference signal [ns] 



2.5 3 3.5 4 4.5 5 
single photoelectron charge [pe] 



Figure 8. Impact of LED calibration: Time (panel a)) and SPE charge (panel b)) distributions before (black 
dashed line) and after (solid red line) applying calibration data. 



Evaluating the calibration runs over three years of data taking, the time offsets of individual 
PMTs are found to change by an average of 0.5 ns from run to run, while the SPE charge calibration 
shows a 10 % spread. The largest values observed are 4 ns in timing and 40 % in SPE charge. For 
about 15 % of the channels, no calibration data is available, as the optical connection via the glass 
fibre is defective. These channels are not corrected for time; their SPE charge is corrected by the 
mean value determined for all PMTs. 

The use of LEDs provides a larger dynamic range of light output intensities compared to the ID 
system that is based on a single laser pulse distributed to the individual PMTs 0]. The OD system 
was designed for a charge calibration of PMTs not only in SPE mode, but also in the region above 
40 photoelectrons (PE), given the much wider dynamic range of PMT pulses expected. However, 
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only time and SPE calibration are applied in regular intervals. 



3. The Muon Veto 

Muons are the only relevant component of cosmic radiation penetrating the rock shielding of about 
3 800m.w.e. over the LNGS underground Halls. At this depth, the muon surface flux is reduced 
to ~1.2 m _2 h _1 . Due to their mean energy of (E^) « 270 GeV flED, most muons travel through the 
detector without any significant loss of energy. 




Figure 9. Definition of muon subsamples. Muons crossing only the Outer Detector are abbreviated as 
ODja's. Muons traversing the Inner Detector (IDji's) are further distinguished in Inner Vessel-crossing muons 
(IVji's) and only Buffer-crossing muons (Bja's). In addition, the position of the auxiliary External Muon 
Tracker (EMT) is shown (fOJ). 

Fig. |9| shows the possible scenarios of muons crossing Borexino. Muons passing through the 
scintillator of the Inner Vessel (IVjli's) create a very large light output as they deposit hundreds of 
MeV of ionization energy. In this case, the limitations of both the PMTs and the electronics cause 
a distortion of the energy information. Still, the visible energy is comparable to point-like events 
of several tens of MeVs. As this energy is far above the 7 Be-v shoulder or even the endpoint of 
the 8 B-v recoil energy spectrum, IVji's are not a serious background for solar neutrino observation. 
The situation changes for muons crossing the ID buffer without touching the IV. Buffer muons 
(Bji's) generate light partly by the Cerenkov effect and partly by the residual scintillation of the 
buffer liquid, despite the DMP quencher dissolved in it [0]. These situations are dangerous, as the 
B|u light output is low enough to mimic a solar neutrino signal. The amount of produced light 
depends almost linearly on the length of the muon trajectory in the buffer. The corresponding 
muon signals cover therefore the whole range of solar neutrino recoil energies. The visible energy 
spectra of muons and point-like scintillation events are presented in fig. [ITJ Here and in the rest 
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of this paper, the energy variable adopted is the number of PMT hits. The relation between this 
variable and the energy actually released in the detector is therefore linear only in the sub-MeV 
range, while at higher energy the effect of multiple photoelectrons being generated in the same 
channel and contributing to the same hit determines deviations from linearity. More details can be 
found in Efl. 

Finally muons crossing the OD only (ODji's) are acquired with a special trigger identification 
flag indicating that no light is present in the ID. Unless otherwise noted these events are not included 
in any analysis. 
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Figure 10. Visible energy spectra (in number of PMT hits) for point-like (dashed red) and muon-like (solid 
blue) events in the Inner Detector, the latter flagged by the OD muon trigger (MTF). 

The probability of a muon to stop in the detector is: 

Z av 

C p = y-Ppc (3-D 

where ppc = 0.88 g/cm 3 is the density of Pseudocumene, A ~ 700m.w.e. is the absorption length 
of cosmic muons for the LNGS rock coverage IfTOl and Z* v is the average path length of the muons in 
the detector. Due to saturation effects following the passage of an IV|i, we are only able to observe 
muon decay of Bji's. In this case Zj^ = 7.2 m and we obtain Pg° p ~ 1 .0%. To check this expectation, 
we select Bja's events that feature two clusters of hits and reject those that can be ascribed to PMT 
after-pulses and other re-clustering effects in the muon signal tail, via the application of specifically 
tuned pulse shape cuts. In the remaining sample the second cluster is therefore due to muon decay. 
This is testified by the profile of the time difference between the two clusters, which correctly shows 
an exponential behavior. The fit returns # 2 /ndf = 158/124 and t = (2. 14 ± 0.05) (Lis. Herein, ndf 
stays for the number of degrees of freedom. An approximate evaluation of the cut efficiencies 
allows us to infer a rate of decaying muons of ^28 cpd, corresponding to ~1% of Bja's, in good 
agreement with the prediction. 

3.1 Muon Identification Criteria 

Muon identification tags concern both the OD and the ID. The OD muon tags search for an increase 
in OD PMT activity above usual dark noise fluctuations; they operate at two different stages of 
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data processing, one hardware, the other software. For the ID we observe that typical muon pulses 
feature longer rise and decay times than neutrino-like scintillation events of comparable energy, 
as - depending on the track length - a muon crossing the buffer needs up to 45 ns to deposit all its 
ionization energy. In addition, the number of detected photons is several orders of magnitude above 
the usual light output of point-like events, prolonging the signal substantially by generating large 
amounts of PMT after pulses. Pulse shape discrimination (PSD) can therefore be used to identify 
muons. 

3.1.1 Muon Trigger Flag (MTF) 

The first tag is a hardware trigger of the OD (see sec. |2.3.4]) , called the Muon Trigger Flag (MTF). 
The trigger condition is given by a significant increase in number of hits in a time span comparable 
to the transit time of a muon. The threshold is set to 6 PMTs firing within a timegate of 150 ns. 

3.1.2 Muon Clustering Flag (MCF) 

The Muon Clustering (MCF) algorithm is an offline software trigger based on the acquired OD 
data. Its trigger condition is a slight refinement of the MTF: two subsets of PMTs are regarded 
separately, distinguishing between hits at PMTs mounted on the SSS and PMTs mounted on the 
floor. The trigger condition is met if 4 PMTs of either subset fired within 150 ns. 

3.1.3 Inner Detector Flag (IDF) 

The Borexino analysis software uses a number of variables that are sensitive to either rise or decay 
time of the acquired pulses. The leading edge of the signal is described by the peak time variable 
t p , the time difference between start and peak of the pulse. The mean time t m is defined as the 
average of the time differences between the start of the pulse and the individual hits in the pulse 
and is sensitive to the overall signal duration. 




visible energy [number of hits] visible energy [number of hits] 

Figure 11. Mean time t m (left panel) and peak time t p (right panel) variables as a function of the visible 
energy of ID events. Blue dots represent events flagged by the OD tag MTF as muon, orange dots correspond 
to point-like events. 

Fig. [TT] shows t m and t p as a function of the visible event energy E v [ s . In the low (100—900 
hits) and medium E Y { S range (900—2 100 hits), t p is the more efficient discrimination parameter. 
The muon discrimination thresholds are chosen to be t p > 40 ns for the low range and t p > 30 ns for 
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the medium one. For high energies, PMT after-pulses begin to contribute significantly to the signal, 
increasing especially t m . In the highest range at E v { s > 2 100 hits, a cut t m > 100 ns is applied. In 
general, the limits were chosen conservatively to avoid the false identification of point-like events 
as muons. The pulse shape Gatti parameter g 9 is used in the low energy region to reject scintillation 
events located in the most upper region of the Inner Buffer (z > 4 m) and arising from a small leak of 
scintillator into this volume. These very peripherical events show t p values larger than for the other 
scintillation events and would therefore be identified as muons, however they are still recognized as 
j8 -like by the Gatti algorithm. In the high energy region the Gatti parameter is instead used to reject 
events due to electronic noise or re-trigger on after-pulses after the actual muon. These events were 
present only in the very first months of data taking, before a hardware protection was introduced, 
and show very high g values that allow an easy identification. 
The criteria for muon identification are summarized in tab. Q. 



Table 2. Criteria for Muon Identification in the Inner Detector. 



Number of Hits 


Mean Time 


Peak Time 
h 


Pulse Shape 
Gatti Parameter g 


100-900 
900-2100 
>2100 


X 
X 

>100ns 


>40ns 
>30ns 

X 


>0.2 (ifz>4m) 

X 

<0.55 



3.2 Muon Tagging Efficiency 

The efficiencies of the available muon tags were investigated using two approaches. First, samples 
of muon events were selected independently of any tagging flag and have been used to find directly 
the identification efficiencies (sec. |3.2.1|) . Second, each flag has been tested against a muon sample 
selected by either of the remaining two flags (sec. |3.2.2l) . 

The analysis of this section focuses on IDji's and includes an energy threshold of 80 hits. 
This value is not higher than the energy threshold applied in any published or foreseen Borexino 
analysis. We exploited all available data of years 2008 and 2009 except runs not passing the 
validation criteria fOl . The remaining data set corresponds to an exposure of ^500d including 
^2-10 7 scintillation events, ^210 6 cosmic muons and ^2.7- 10 4 events arising from the CNGS 
beam (see app. 

In all these studies, the selected test samples contained very small (if any) impurities due to 
non-muon events. When sizable, such impurities have been subtracted from the sample; however, 
we conservatively regard the resulting efficiencies as lower limits. 

All errors reported on efficiencies throughout the section have been computed at 95% C.L. 
using Clopper Pearson Intervals. 

9 The Gatti parameter [Oil is computed from a pulse shape analysis. Its value is deduced from an algorithm imple- 
menting an optimum filter technique and developed for the discrimination of a-like (g > 0) and jS-like (g < 0) events. 
Reference shapes originating from MC studies in combination with intrinsic a and /3 contaminants of the scintillator 
such as 214 Bi and 214 Po whose fast coincidence allows a very clean sample selection [Pi . 
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3.2.1 Direct Efficiency Determination 

The key requirement in selecting a test muon sample is that the flag being tested is not dependent 
on parameters used for the selection. As detailed below, high-energetic ID events present in normal 
neutrino runs provide an appropriate test sample for the OD related flags, while events on time with 
CNGS Vju beam spills provide a reliable test sample for IDF. 

Direct OD flags efficiency via high energy ID events 

Fig. [T^ shows the ID energy spectra of events both flagged and not flagged by the OD MTF tag. 
As can be seen, at high visible energies, scintillation events in the ID are predominantly caused 
by IVjLi's. We selected events with more than 7 000 hits, corresponding to a visible energy of 
^200 MeV or ^10 5 photoelectrons. This estimate can be derived from the fraction of after-pulses 
in the signal (~5 000) and from the average probability of after-pulsing (5 %). This energy is higher 
than the most energetic solar neutrinos (the hep-v's with an end-point energy of 18.8 MeV), and far 
above the maximum energies of the radioactive decays in the natural U/Th chains ( 208 T1 features a 
Q-value of 5.0 MeV). 

The reason why we regard this choice as valid is that the efficiencies of the OD flags are 
independent of the visible energy detected in the ID: the PMT coverage of the OD shows no relevant 
gaps to the cosmic muon flux; all such muons — independently of the exact track geometry — 
have to cross several meters of water before entering the steel sphere, thereby creating comparable 
amounts of Cerenkov light. The result of this study can therefore be assumed to hold over the 
whole energy range of ID events. 

Impurities in the selected sample could be caused by fast neutrons or by ID-contained atmo- 
spheric neutrino events. Although we cannot measure this contamination of the sample we expect 
such events to account for a very small fraction and we neglect them. 

The resulting efficiencies are £ M tf > (99.25 ±0.02) % for MTF and £ M cf > (99.28 ±0.02) % 
for MCF (see first column of the "Direct Efficiencies" section in tab. 0). 

Direct ID flag efficiency via CNGS muon events 

The CNGS neutrino beam (see details in app. 0) provides a sample of muon events created by 
Vjj charged current reactions in the upstream rock or the detector. As these events can be easily 
assigned to the beam spills via their GPS time stamp, they are in principle an ideal test sample for 
the muon tagging flags. However, the OD is optimized for vetoing and tracking the cosmic muons 
which arrive preferentially around the vertical direction. On the contrary, the CNGS neutrino 
beam arrives at the detector at an angle of about 3.2° below the horizon, and therefore generates 
mostly horizontal muon tracks. As a result, the distribution of the PMTs is not ideal (especially 
for muons crossing the non-instrumented lower quarter of the sphere) which leads to a decreased 
veto efficiency. The ID instead, due to its spherical symmetry, features an ability for muon flagging 
which is independent of the track direction. Upon these considerations, CNGS muon events are not 
a good test sample for OD-related flags, while they are perfectly valid for IDF. 

A contamination of the sample arises from 14 C and 210 Po decay events which can occur in 
accidental coincidence with the beam spills. In app. |A| we have computed that < 0.06% of the 
sample events considered here can be regarded as accidental coincidences. Although it is a very 
small effect, these events have been subtracted from the test sample. 
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We have tested the IDF tag in four energy regions. The results are shown in the second column 
of the "Direct Efficiencies" section of tab. |3|. The general increase of efficiency with E v { s is expected 
as the larger number of registered PMT hits improves the temporal pulse shape resolution. The 
global energy-independent efficiency was determined to be > 95.86^1^%. 

Table 3. Efficiencies of the available muon identification tags. All numbers are lower limits as the presence 
of non muon events in the test samples cannot be excluded. Errors not shown if smaller than 10 -4 . 



Tag 
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Direct Efficiencies 


Mutual Efficiencies 
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Figure 12. IDF tagging efficiency against MTF- selected and against CNGS -selected events as a function 
of visible energy. 



3.2.2 Mutual Efficiency Determination 

The efficiencies of the three tags can also be derived by mutual comparison. Each of the three 
flags has been tested against muon samples selected using either of the two remaining flags. The 
results of the mutual tests of MTF and MCF to each other are included for completeness and 
should be handled with care as the two flags share the same set of PMTs and are therefore not truly 
uncorrected. However, tests seeing IDF on the one side and MTF/MCF on the other are fully valid 
as the two detectors are completely independent. 
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Impurities in the samples arise here from point-like events due to potential "over-efficiency" 
of the reference flag. MTF and MCF are in principle free of over-efficiencies as there is no light 
connection between ID and OD. However, small light leaks in the external steel dome increase the 
rate of OD|i-like events up to ~2.5- 10 4 cpd and this can cause accidental coincidences with point- 
like ID events. The number of such events has been statistically evaluated; we have found that for 
MCF (MTF) the "over-efficiency" above the 80 hits threshold are as low as ~10" 5 (-10" 6 ). The 
over-efficiency of IDF is not known precisely, however from the results of this section it is limited 
to below 0.7%. 

The resulting efficiencies are shown in section "Mutual Efficiencies" of tab. 0. The tests versus 
IDF imply efficiencies > (99.33 ± 0.01) % for MTF and MCF. IDF on the other hand shows a 
good performance, > (98.90 ±0.01) %, if integrated over the whole muon spectrum. However, in 
agreement with the results of the direct efficiency study, discrimination power decreases for low 
visible energies (fig. |T^|) . 

3.2.3 Overall Veto Efficiency 

We have evaluated efficiency of the available muon tagging methods both directly and from mutual 
comparison obtaining consistent results that are summarized in tab. |3|. 

In all Borexino analyses, muons are suppressed from the data set by the application of both 
OD- and ID-related flags, MCF being the preferred choice when available to combine with IDF. 
MTF is used in the small fraction of data when MCF is not available due to DAQ inefficiencies. 
The muon veto efficiency from the application of the combined tags is at least 99.992 %. 

4. Muon Track Reconstruction 

The primary goal of muon tracking in Borexino is to enhance the veto of cosmogenic radioisotopes 
by the introduction of geometrical cuts. Muons produce unstable nuclei by spallation processes 
on their trajectory through the scintillator and the decay of such nuclei can mimic the expected 
signal. For short-lived radioisotopes, it is usually sufficient to veto all events in the wake of a 
muon for a duration corresponding to several half-lifes of the given isotope. This kind of veto 
is sufficient for 7 Be-v and geo-neutrino analyses because the cosmogenic background in these 
cases is due to nuclei whose half-life is of the order of a few seconds or less. However, when 
the half-life of the dangerous nuclei is longer, the induced dead time can become very large. To 
overcome this problem, the veto of events can be applied only to a cylindrical region around the 
muon track instead of the whole detection volume. While information on the muon track was 
already exploited in the 8 B-v analyses [0], it is of uttermost importance for the CNO and pep- 
V search: here, the most prominent background arises from cosmogenic 11 C that surpasses the 
expected neutrino rate by an order of magnitude. n C features a half-life of 20.38 min, so a complete 
veto of the detector for several half-lifes following each muon would result in a dead exposure close 
to 100 %. Fortunately, 11 C production is in 95 % of all cases accompanied by a knock-out neutron 
from 12 C. The n C background events can be identified by a threefold coincidence relying on spatial 
and time information of the muon, the neutron (see sec. |5|), and the decay signal of the radionuclide 
as described in more detail in flFUl . 
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Moreover, muon tracking allows an accurate investigation of both cosmogenic neutron back- 
ground and other muon-induced production of radioisotopes occurring in Borexino. Such a study 
is in preparation and is of particular value not only for the Borexino analyses, but also for other low 
background experiments located at LNGS or at other underground laboratories of similar depth. 

4.1 Reconstruction Algorithms 

We have developed two independent tracking algorithms. The OD tracking is based on the entry 
and exit points generated by the Cerenkov cone of muons in water. The ID tracking relies on the 
distinctive photon arrival time pattern created by the scintillation light front emitted by the muon 
track. The results of both ID and OD tracking algorithms are later combined in a best fit global 
tracking. 

4.1.1 Outer Detector Tracking 

The algorithm is based on the identification of the disc-like activation profile of OD PMTs gener- 
ated by IDji's: the Cerenkov cone around the track is registered by the PMTs close to the muon 
entry point EP (usually on the SSS) and with some time delay by PMTs near the exit point XP (ei- 
ther on SSS or floor of the water tank, depending on track inclination). High-energy cosmic muons 
are highly relativistic. As the Cerenkov light is propagating at the speed of light in water (0.75c), 
the muon arrives at the PMT-equipped surface first, followed by its Cerenkov light cone. If the 
track is perpendicular to the surface, the first light is detected by the PMTs nearby the penetration 
point. As light is emitted at a fixed Cerenkov angle a c = 41°, the time profile of PMT hits follows 
circles or (depending on the track inclination) ellipses around the penetration point, from early hits 
in the center to the late hits at the periphery. The number of photons per PMT decreases with the 
distance to the penetration point. Fig. |T3| shows a view of the typical PMT activation profile. 
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Figure 13. Activation profile created by the Cerenkov cone of a muon that is crossing the curved surface 
of the SSS: the size of the dots indicates the number of photons registered by a PMT, while colors code the 
photon arrival time (y Metered) . 

The tracking algorithm uses clusters of hits associated in time and space to identify EP and 
XP. PMT hits on the SSS and on the floor are treated separately. The spatial distance between two 
hits on the SSS (floor) has to be Sr < 3 m (4 m) in order to affiliate the pair to the same cluster. 
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This condition corresponds to about twice the mean distance between neighboring PMTs, which 
is 1.7 m (2.2 m). In addition, the time difference between these hits has to be smaller than 20 ns, 
which is generous compared to the expected time difference between neighboring PMTs of 4.9 ns 
(6.9 ns). In this way, separate clustering of light reflected back at the original penetration point is 
suppressed. 

The charge bary centers R\> c of clusters are computed to serve as future entry and exit points. 
The coordinates f{ of the affiliated hits are weighted by Wi considering both the signal charge qt and 
the PMT hit time t t : 

4c = with wi(q u ti) = qt • e-(«-*>/ T (4.1) 

with t = 20 ns and to the time of the first hit affiliated to the cluster. The barycenters are calculated 
in Euclidean coordinates, but radially shifted to the radius of the OD PMT photocathodes to take 
into account the SSS surface curvature. The deviation of the weighted mean is used as statistical 
uncertainty. 

For an ideal track, only two clusters caused by EP and XP belonging to a single track are 
present, making no further selection necessary. However, light reflected by the Tyvek foils creates 
in most cases additional fake clusters that have to be rejected. A large fraction of these clusters 
can be discriminated by setting a minimum for the integral number of photoelectrons as most of 
the light intensity is concentrated in the clusters containing direct light. The main criteria applied 
to identify the OD entry cluster are its total number of photoelectrons and its synchronization with 
the start time of the corresponding ID event. Based on the identified EP, the exit cluster is found 
by checking the time-of-flight difference that corresponds to the trajectory of the muon traveling at 
speed of light in vacuum. Moreover, a reconstructed track is only considered valid if it crossed the 
sphere, as due to the synchronous trigger of ID and OD only IDji's are considered. 

There are events where only the EP can be identified: if the muon stops inside the ID, there is 
no XP. Alternatively, if the track leaves the ID through the optically inactive lower hemisphere at 
a flat angle, it produces only a small amount of light that is sometimes too faint to be detected. In 
these cases, it is vital to avoid the identification of a reflection light cluster as XP. 

4.1.2 Inner Detector Tracking 

The ID muon reconstruction is based on several successive fits to the distinct photon arrival time 
patterns generated by the tracks. The muon entry point EP is found by fitting the arrival time 
distribution of the first hits detected inside the ID. We then introduce the plane P exo , defined by 
the muon entry and exit points as well as the SSS origin. P exo is a symmetry plane that divides 
the photon arrival pattern of the PMTs in two identical halves. It is identified by several fits to the 
arrival time patterns along rings of PMTs in rising distance from EP. As the muon track lies on P exo , 
both EP and the exit point XP are in this plane (fig. [P]). Finally, a fit to the time distribution of the 
PMTs in P eX o is used to find the discontinuity in the time pattern around XP. This local minimum 
at XP is generated by the tip of the muon light cone penetrating the sphere surface. 

In fact, the light front generated by a muon in the ID is described by a superposition of scin- 
tillation and Cerenkov light. While scintillation light is clearly dominant for the part of a track 
located in the IV, the DMP added to the buffer suppresses the scintillation light output of the PC 
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Figure 14. Geometry defining the muon track: the coordinate system was rotated to shift the entry point 
EP to the zenith of the SSS. The symmetry plane P exo of the photon arrival time distribution is defined by the 
sphere origin O, EP and the exit point XP. The azimuthal coordinate (p x of XP is determined by fits to the 
arrival time distributions on horizontal circles of PMTs (parallel to the red dashed circle). The zenith angle 
is found by a fit to the arrival times along the intersection of P exo and the SSS surface (see fig. |17| ). 

JJ. For Bu's and for the buffer portion of the track of IVu's, light is generated in comparable 
amounts by scintillation and the Cerenkov effect. However, the resulting shape of the light front 
will be similar in IV and buffer: as depicted in fig. [T3], the superposition of spherical waves of scin- 
tillation light along the muon track forms a conical light front centered around the muon trajectory. 
The opening angle is identical to the Cerenkov angle as it depends merely on the light speed in the 
scintillator medium. In addition to the cone, a backward-running spherical light front is caused by 
the scintillation light only. 

Entry Point Reconstruction 

The backward-running light reaches the SSS first at the muon entry point EP. For a first approxima- 
tion of the EP position, the ID tracking analyzes the time pattern of the hits detected in the first 5 ns 
of the muon event, corresponding to the first 70 cm of the track inside the ID: similar to the OD, the 
quenched scintillation light of the buffer produces ellipses of isochronous PMT hits for which EP 
lies approximately in a focus (fig. |T^). As the charge information of ID PMTs is not reliable due to 
the large muon light output, only timing is used to compute the barycenter according to eq. [O]. 

Based on the first guess, the time distribution of hits around EP is fitted with a function t e ((j>) 
optimized to find the minimum arrival time that corresponds to EP: for this, the first-approximation 
barycenter is shifted by rotation to the equatorial plane (azimuth angle <p = 0, zenith angle & = n/2) 
in order to avoid the topological effects of the spherical coordinates that arise mainly on the poles. 
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Figure 15. A muon traversing the ID creates Figure 16. A muon entering the ID creates ellipses of 
light by scintillation processes. The superposi- isochronous PMT hits. The entry point is in one focus, 
tion of spherical waves along the track creates The very first hits (symbolized by the violet area) should 
a light front that is in forward direction identi- be almost centered on the entry point position. The relative 
cal to a regular Cerenkov cone. In addition, a offset of the bary centers for later isochronous hits bears in- 
spherical light front propagating backwards is formation on the inclination vector of the muon track d. It 
created. is used in an alternative algorithm described in IfTHl . 

The empirical fit function 

I a + • sin(c + (<p' - <p' e )) else 

is of V-like shape. It allows for different curvatures and slopes a ± of the arms in order to 
accommodate the curved surface geometry and a possible relative inclination of the track that 
results in differences in the light propagation time. Also a small time offset a is included. Fig. [17 



shows the time distribution of hits dependent on the angle relative to the barycenter. Plots a) and 
b) correspond to two perpendicular orientations along the sphere surface that correspond to the <p' 
and planes in the transformed coordinates. The deviations <p' e (and ^) from <p' = = 0) that 
can be extracted from the fits are transformed back to regular Borexino coordinates to obtain the 
corrected position of EP. 

Finding the Symmetry Plane 

To identify the muon exit point XP, the algorithm searches for the symmetry plane P exo defined 
above. To find it, the PMT coordinates are again rotated, this time to shift EP to the zenith. The 
corresponding geometry is shown in fig. 0. If in these coordinates a horizontal plane of PMTs 
is chosen and the arrival time is plotted as a function of the azimuth of the transformed PMT 
positions, the observed time pattern is again of V-like shape. The minimum is located at the angle 
<p' x corresponding to P exo in the transformed coordinates. The fit function 



t x (q>') = c'l • A /4+^ I 2 -2^s^i(cos(p / cos^ + sin9 / sin(p;)+/o (4.3) 

reflects the symmetric dependence of the arrival times on (p f . It is inspired by the time distribution 
expected if the track crossed the chosen plane of PMTs exactly perpendicularly: in this case, Rs 
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Figure 17. The arrival time pattern at the entry point registered by ID PMTs as a function of zenith (a) 
and azimuth (p f (b) angles corresponding to PMT positions on a geodetic circle around the SSS. The two 
circles intersect at the supposed positions of the entry point at (p' = = 0. 



and R\ would denote the radius of the sphere and the impact parameter of the muon track relative 
to the sphere center, respectively, while c pc denotes the light speed in liquid scintillator. 

The algorithm performs this fit on twelve horizontal planes of PMTs in rising distance below 
EP, each returning a best-fit value for (p x . Finally, the mean of these results is used, weighting each 
result with the fit uncertainty. In principle, this algorithm returns a very precise result for (p x as the 
information of almost all of the PMTs is taken into account. However, if EP position is off by a 
certain angle relative to P exo , this will cause a sine-like distortion of the (p x fit results (see below). 

Exit Point Identification 

To find the second coordinate angle the PMT circle defined by P exo is selected. The correspond- 
ing arrival times are drawn in fig. p~8| a as a function of the (transformed) zenith angle . Here, the 
EP is located at ^ = 0. The arrival time distribution shows a discontinuity in the rise for E [0, n] . 
The kink is caused by the projection of the Cerenkov-like light cone on the surface defined by the 
SSS, the local minimum corresponding to The fit function 



t x (#) = i 5 + *o (4.4) 

cyjl - £ • (sin(arccos(i) - - ^|) +nsin(±|tf' - ^|)) 

is derived from an analytical time of flight calculation based on the sphere radius Rs, the speed 
of light c and the refractive index n. The function assumes light emission from the muon track 
under an angle corresponding to the Cerenkov condition and provides an excellent fit. However, 
the function is only valid close to the true value of for other regions, the observed arrival 
times are changed due to the geometrical cut-off caused by the sphere. Therefore, the fit routine is 
performed only on a limited interval surrounding a first prediction of $ x based on the visible light 
output: fig. [T9| shows the basic relation between the number of hits registered in the ID and the track 
impact parameter R\ determined with the OD tracking algorithm. The fit routine is iterated several 
times, each time starting from the $ x computed in the previous step and progressively narrowing 
the fit range at each step. 
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This method works best above 4000 hits, which corresponds to tracks with an impact param- 
eter Ri of less than 4.5 m (fig. [L9|). In the region from 2500 to 4000 hits (Ri e [4.5 m; 6m]), the 
scatter in the arrival times becomes larger so that the distribution has to be averaged over angular 
bins of A# = 7T/30 in width. In addition, a slightly modified fit function that stresses the disconti- 
nuity around XP has to be applied in order to obtain reliable results. Below 2500 hits, $' e and id-' x 
lie too close to be discerned in the fit, and id-' x collapses to 0. In this case, only EP coordinates are 
later on used in the global track (see below). 




d" [rad] 9* [rad] 

Figure 18. Photon arrival time distribution around the exit point XP: a) distribution for PMTs on the 
symmetry plane of the sphere P exo (see fig. [Hf); b) distribution (p f for PMTs corresponding to a geodetic 
circle intersecting at 
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Figure 19. The impact parameter determined by OD tracking versus the number of ID hits in the muon 
event. Two main bands are visible in the regions below and above 7000 hits: they are caused by Bja's 
and IVja's respectively. The IV band is strongly bent due to the relation of path length and of the impact 
parameter. 



Finally, the fit on t x {cp') is repeated based on the effective function described by eq. |4~2| . The 
result is shown in fig. [T8|b. This last step is used to apply a small correction to <p' x in cases in 



which the symmetry plane P exo determined before is inclined relative to the true muon track. XP 
coordinates are then obtained by transforming back (<jp£, to the original reference frame. 

4.1.3 Global Tracking 

Combining the information provided by OD and ID muon tracking, a global track can be defined 
representing the best fit to the (up to four) data points available for a single track. Different from 
the other tracks, the information on the global track is not saved in the form of points on the track 
but using a set of four fit parameters. The three dimensional track T{x) is parameterized as : 
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The reduced notation is y(x) = a + j8x, z(x) = y+ 8x. The latter equations already indicate that the 
3D fit can be split into two independent 2D fits to projections of the track on the xy and xz coor- 
dinate planes. To preserve the information on the direction of motion of the muon, a fifth variable 
indicating the track orientation is necessary. Fig. \ZQ shows the resulting track in a graphical repre- 
sentation. The four track points including their uncertainties (indicated by boxes) are shown along 
with the fitted track; in addition, time and charge profiles provided by ID and OD are included. 




Figure 20. Graphical description of the global muon tracking. The track points provided by ID and 
OD reconstruction algorithms are shown, the boxes indicating their uncertainties. Also included is the 
information on the ID PMT hit time distribution (small colored dots) and the OD hits (colored circles). 

There are a number of cases in which 3 out of 4 points are reasonably aligned while a fourth 
point seems to be misreconstructed. If the reduced % 2 value is # 2 /ndf > 3, four additional fits are 
performed that each omit one of the track points and the fit with the lowest # 2 /ndf is chosen. 
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4.2 Tracking Performance 

The parameters entering the tracking algorithms were tuned using a MC simulation of muon events. 
The evaluation of the accuracy of the tracking instead has been performed on real data. Several as- 
pects are evaluated: 

The efficiency of the tracking can be tested by the ratio of cosmic muon events that feature re- 
constructed tracks: complete tracks with both entry and exit point are found for 98.5 % of all muon 
events for the OD tracking. The efficiency of the ID tracking is lower, 87.9 %. The reduction is 
mostly due to the minimum energy requirement of 2500 hits; above this threshold, the efficiency 
is 96.9 %. The global tracking algorithm returns a fitted track in 96.4 % of all muon events, cor- 
responding to 99.3 % of the cases in which at least two points were found by ID or OD tracking. 
These values on their own have no further implication concerning the quality of the reconstruction. 

The angular resolution of the tracking can be determined by utilizing CNGS muon events, gen- 
erated in the rock upstream of the Borexino detector and crossing both Borexino and OPERA 
DH 0,01]. A track-by-track comparison between the two experiments is possible, returning reso- 
lutions of 3° for OD/global tracking and 5° for ID tracking. In addition, a high-resolution External 
Muon Tracker (EMT) has been mounted on top of the Borexino steel dome in late 2009 (see fig. ^), 
providing a small sample of well-determined muon tracks crossing Borexino. Statistics of this 
analysis is still limited, but the resulting resolution is comparable to the one obtained from CNGS 
muons. The results of CNGS and EMT analyses are complementary as they test horizontal and 
vertical muon tracks, respectively. 

The lateral resolution (perpendicular to the muon track) can be derived from a track-by-track com- 
parison with OPERA, returning a resolution of 40-50 cm (35 cm for IVji's). Statistics for EMT are 
in this case too low to return meaningful results. An alternative is the analysis of the distribution 
of cosmogenic neutrons relative to their parent muon track. Here, a resolution of ^35 cm for IVji's 
is obtained. 

In the following, CNGS, EMT and cosmic results are discussed in more detail. While the perfor- 
mances of the individual tracking algorithms depend in detail on the applied test, consistent values 
can be found. Generally speaking, vertical and through-going muon tracks are best reconstructed 
by the OD and global trackings, while horizontal and (semi-)contained tracks are reproduced best 
by the ID algorithm. 

4.2.1 CNGS Muons 

As described in app. the CNGS beam arrives at the LNGS at an inclination angle of 6 = 93.2°. 
Relative to CNGS, CERN is located at an astronomical azimuth of <p = 305°, almost exactly coin- 
ciding with the orientation of Hall C. Fig. shows the azimuth and inclination angle distributions 
obtained by Borexino (and OPERA) tracking. 

A test of the tracking angular resolution is possible by an evaluation of the CNGS events that 
were reconstructed both in Borexino and in OPERA. The angular resolution for muon tracks in 



-25- 



- OPERA track 

- global track 

■ OD track 

■ ID track 



50 



100 



150 



200 



250 




140 160 180 
inclination 9 [°] 



Figure 21. Angular distributions of muons produced by the CNGS beam. Tracks in Borexino are 
required to cross both ID and OD. The distributions reconstructed by OPERA are shown for comparison. 



OPERA is of the order of mrad and by far superior to the performance of the Borexino tracking 
IfTgll. Therefore, the intermediate angle a between the reconstructed track orientation vectors of 
both experiments corresponds almost exactly to the angular resolution of the Borexino tracks: 

cos a = /i((p, 6) • v(0,#) = sin0 sin#(cos<pcos0 + sin<psin0) + cos0cos# (4.6) 

where /i(<p, 0) and v(0, #) represent the track orientation vectors. 

In the period from 2008 to 2010, 2797 events featuring tracks in both Borexino and OPERA 
detectors are found (^5 % of the overall number of Borexino ID on-time events). Several selec- 
tion criteria are applied to allow a comparison: in Borexino, events are required to produce light 
synchronously in ID and OD, in order to obtain a sample of through-going muon tracks, generated 
in the rock upstream of the detector. In OPERA, only muon tracks featuring a momentum greater 
than 20 GeV are selected. This reduces the effect of muon forward-scattering that might otherwise 
bend some of the muon tracks on their passage from Borexino and OPERA. The resulting samples 
of coincident muons feature a few hundred events. 

Fig. ^ shows the distribution of a as defined in eq. £|15| for the comparison of OD and OPERA 
tracking. Each entry is weighted with 1 / sin a to take the spherical geometry into account. The 
distribution can be fitted by the sum of a Gaussian centered on a = and a constant term that 
accounts for misaligned tracks for larger values of a. The latter are mostly caused by muons 
emerging from charged-current interactions inside the Borexino ID and missing a reliable entry 
point for both ID and OD tracking. 

Tab. £| lists the fit results for the Borexino tracking algorithms. The angular resolutions deter- 
mined in this way range from 3° (OD, global) to 5° (ID tracking). If only IVjli's are selected, the 
resolution improves to ^2.5°. 

As not only the track direction but also the absolute track coordinates are known for both 
Borexino and OPERA, a determination of the lateral resolution is possible. For this, the distribution 
of relative distances of the track penetration points in the Borexino yz-plane (perpendicular to the 
CNGS beam direction) can be used. Fig. shows the y-distance distribution between Borexino 
global and OPERA tracks. The center was artificially set to 0, as the relative coordinates of the 
two detectors are known only at an accuracy of 1 m. However, the distance distribution can be 
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X 2 /ndf 22.7/24 

norm 209+ 14.6 

mean 0.024510.0315 

a 0.437+0.028 



X 2 /ndf 11.7/12 

norm 52.6 + 7.3 

mean 0.0849 + 0.0478 

a 0.325 ± 0.046 
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Figure 22. Intermediate angle a between the direc- 
tion vectors of muon tracks reconstructed in OPERA 
and Borexino OD tracking. The width of the Gaus- 
sian fit corresponds to the OD angular resolution. 



Figure 23. y-distance distribution between recon- 
structed tracks in the Borexino yz-plane, compar- 
ing OPERA and Borexino global tracking. The lat- 
eral resolution a y improves if only IVu's are selected 
(shaded red area). 



approximated by a Gaussian, the width o y corresponding to the lateral resolution. For the vertical 
z direction, the same fit can be applied, returning o z . 

The resolution values for all trackings are reported in tab. [|. Considering all IDji's, there is 
no difference in a y and a z for ID tracking, but large values for a z are observed for the global and 
especially for OD tracking. This effect is caused by the lower non-instrumented section of the SSS 
outer surface - horizontal OD tracks passing this section will be reconstructed either too close to the 
equator or to the OD floor. Obviously, also the quality of the global tracking is affected. However, 
the effect vanishes if only tracks passing the PMT-equipped upper part of the SSS are selected. If 
only IVji's are selected (more than 7000 hits in Borexino), the lateral resolution improves in both y 
and z to ^35 cm for all trackings, a y « a z w Oi at . 



Table 4. Fit results for angular and lateral resolution obtained from CNGS muon tracks observed by 
Borexino and OPERA. 



IDjj's 


<T«[°] 


o y [cm] 


o z [cm] 


OD tracking 


3.51±0.13 


53±4 


198±13 


ID tracking 


5.13±0.25 


42±3 


40±5 


GL tracking 


2.83±0.13 


44±3 


87±12 


IVu's 


<T«[°] 


o y [cm] 


o z [cm] 


OD tracking 


3.01±0.15 


28±7 


45 ±7 


ID tracking 


2.44±0.19 


36±5 


31±6 


global tracking 


2.31±0.13 


33±5 


35±7 



4.2.2 External Muon Tracker 

The external muon tracker (EMT) is a support system allowing the determination of the coordinates 
of a subset of muon tracks at high precision in order to test the Borexino reconstruction algorithms. 
Originally developed as a prototype of the OPERA muon tracker, the EMT is based on layers of 
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drift tubes. Whereas in the original setup only a 2D track reconstruction was possible, the tracker 
has been modified for a three dimensional tracking by crossing the modules. It has been mounted on 
a platform above the Borexino steel dome. The location was chosen 4.6 m off the vertical symmetry 
axis of the detector and above the dome in order to sample muons crossing all sub- volumes of the 
detector (fig. g). 
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Figure 24. Sample event from the EMT track display. Open circles correspond to a head-on view of the drift 
tubes in upper (x) and lower (y) layer. The filled colored circles indicate the observed drift times, from which 
the track projections in the yz- and xz-plane are derived. The overlay shows a schematic three-dimensional 
view of the spatial arrangement of the EMT drift tubes. 



As depicted in fig. 0, the EMT modules are arranged in two arrays of drift tubes that are 
oriented perpendicular to each other, one above the other with an interspace of a few centimeters 
caused by the endplates of the modules. Each array contains 96 closely packed drift tubes stacked 
in four layers of 24 tubes. Each tube is 1.05 m in length and 3.8 cm in diameter, and is equipped 
with a gold plated tungsten wire of 45 |Lim diameter in the center as an anode. The total area covered 
is close to 1 m 2 . This corresponds to an expected muon rate of ^25 events per day. Considering the 
position of the EMT with respect to Borexino, the cosmic muon flux and angular distribution in hall 
C, we estimate the daily number of muons crossing both detectors: An average of 9.6 muons/day 
is expected to cross both the EMT and the OD, of which 6.1 also cross the ID and 2.5 the IV. 

The drift tubes are operated with an 80:20-mixture of argon and carbon dioxide as a drift gas. 
Spatial reconstruction relies on the drift times of electrons caused by the muon through ionization 
of the drift gas. The drift time in a hit tube depends on the radial distance of the muon track from 
the axial wire. Each array of tubes provides 2D information on the muon track. Combined, the 
system provides a three-dimensional track. 

The signals of all tubes are read out by a customized TDC board which is also used in OPERA, 
connected to a PC close to the EMT. Data are written to its local hard drive. The trigger is formed 
by two layers of plastic scintillator panels of 2 cm thickness that are read out by 8 PMTs. The 
trigger is issued in case of a coincidence signal in both planes. The trigger logic is based on NIM 
modules. Low and high voltages are also generated on-site in separate modules. A copy of the 
logic trigger signal is sent to the Borexino main DAQ via a 30 m long coaxial cable. This allows us 
to relate corresponding events of the EMT and Borexino by time synchronization. 

Laboratory tests give an angular resolution better than 5mrad and a spatial resolution of 
^300 jam, by far surpassing the resolution achieved for the Borexino tracking algorithms. The 
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position relative to the detector center was determined with an accuracy of ^1 cm. The EMT was 
installed in late 2009, and has been taking data since spring 2010. However, due to the relatively 
low efficiency of the current trigger system the available statistics is rather low. A replacement of 
the trigger system with a higher efficiency one is foreseen in early 201 1. 

As in the case of CNGS muons, the angular resolution a a of the tracking algorithms is deter- 
mined by a track-by-track comparison, relying on the intermediate angle a of eq. fO|). A typical 
histogram for the ID tracking is shown in fig. ^5], the fit results are listed tab. |5|. In spite of the low 
statistics, the obtained values for a a are comparable to the CNGS results presented above. 



X 2 /ndf 12.36/14 

norm 3066 ± 564.7 

a 4.528 + 0.5699 

const 2.067 + 0.9272 
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Figure 25. Intermediate angle a between the direction vec- 
tors of muon tracks reconstructed by the EMT and Borexino ID 
tracking. The width of the Gaussian fit corresponds to the ID 
angular resolution. 



Table 5. Angular resolution as obtained 
by the comparison of Borexino and EMT 
trackings. N is the number of tracks avail- 
able for comparison. 





O a [°] N 


OD tracking 
ID tracking 
GL tracking 


3.9±0.3 125 
4.5±0.6 82 
3.4±0.1 125 



4.2.3 Cosmic Muons 

Cosmic muons are by far the largest muon sample available in Borexino. They allow a qualitative 
test of the angular resolution of the tracking by comparison of the obtained angular distribution to 
the earlier results of the MACRO detector lfT911. Moreover, cosmic muon tracks reconstructed as 
upward-going make it possible to determine the fraction of misreconstructed events. 

The topology of the Gran Sasso mountain massive is imprinted on the cosmic muon flux 
detected at LNGS. The underground halls were dug right under the summit of Monte Aquila, 
which is part of a small mountain ridge on the northern verge of Campo Imperatore, extending 
from West to East. Unlike in deep underground mines in which most of the muons arrive from 
straight above (cos 0=1, being the zenith angle), the reduced shielding towards the mountain 
flanks causes additional contributions at lower values of cos0. The mountain profile (see app. A 
of (USD) reflects in a complicated angular distribution. Fig. shows the measured distributions for 
azimuth (p and inclination angle for cosmic muon tracks detected by Borexino. For comparison, 
the angular distributions measured by the MACRO experiment IfTPIl are indicated as grey-shaded 
areas. 

Only muon tracks reconstructed at cos0 > 0.5 were selected for the shown azimuth distri- 
bution. This corresponds to the limited inclination angle acceptance of MACRO and eases the 
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comparison. For both MACRO and Borexino the tracks peak around (p=30° (N) and <p=170° (S). 
The cos 9 distribution shows a widening in the Borexino trackings compared to MACRO. This 
indicates an inferior track reconstruction in Borexino, OD and global tracking performing slightly 
better than ID tracking. 




azimuth angle cp [°] inclination cose 

Figure 26. Azimuth (p and inclination 6 distributions for cosmic muon tracks detected at the LNGS. The 
colored curves indicate the result of the three Borexino tracking modules. The result of MACRO lfT91l is the 
grey shaded area. 

All three tracking algorithms show a low but relatively constant number of tracks at cos 0<O. 
Upward-going muons are supposed to be very rare as they are solely produced by antipodal at- 
mospheric neutrinos. The tracks at cos 0<O are far too numerous and must be regarded to a large 
extent to be mis-reconstructed. Tab. § shows the fraction of upward-going tracks in dependence 
of the visible energy. Overall, about 2-5 % of tracks are reconstructed as rising. They are mostly 
caused by a misalignment of EPs and XPs for Bjli's crossing the SSS close to its edge, as in this case 
both spacial and temporal alignment of the penetration points is difficult. Therefore, the upward- 
going contribution is reduced to 1-2 % if only IVju's are regarded. 

Table 6. Percentage of mis-reconstructed (upward-going) muons for all IDji's and IVji's only. 





IDjli's IV|i's 


OD tracking 
ID tracking 
GL tracking 


2.0% 1.2% 
5.2% 0.8% 
3.2% 2.1% 



4.2.4 Cosmogenic Neutrons 

Cosmogenic neutrons are visible in Borexino via the y-rays produced after their capture as ex- 
plained in more details in sec. [5]. The analysis of the distance x of cosmogenic neutrons to their 
parent muon tracks is interesting for two reasons. First, the distance distribution of a large neutron 
sample can be used to determine the mean free path of a cosmic neutron in the scintillator. This 
study is in preparation and will be presented as a separate paper. Second, for short distances, the 
x distribution provides information on the precision of the track reconstruction. We focus here on 
this second outcome of the analysis. 
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Figure 27. Distribution of the distance x of neutrons to their parent muons as derived from the global 
tracking. Crosses represent data points, the line corresponds to the MC result. 

Fig. ^7| shows the x distribution derived from global tracking: the neutron sample was selected 
requiring a minimum visible energy of 700 hits, corresponding to neutron capture events on hy- 
drogen inside the FV (fig. ^). This minimizes the effect of spatial uncertainty introduced by the 
reconstruction of the neutron capture vertex but also distorts the information on the mean free path 
of the neutrons by selecting only non- saturating (i. e. non- showering) muons, as explained in more 
detail in sec. I53I . 

This distribution can be reproduced by a Monte Carlo simulation that takes into account the 
muon lateral resolution cfy, the neutron mean free path A, and the uncertainty in the reconstruction 
of the light barycenter of the 2.2 MeV y-rays (o n ). The muon lateral resolution function is com- 
posed of a Gaussian uncertainty a and a constant asymmetric offset ji that can be associated with 
radial shifts of the reconstructed tracks. For the neutron, a n is dominated by the reconstruction 
uncertainty and mean free path of the y-ray (<7 7 « 23 cm was obtained by simulations), but the 
diffusion of the thermalized neutron before capture Odiff ~ 4 cm is also considered lEDl. The best 
agreement to measured Borexino data is obtained by varying the parameters <7, /I, A and the overall 
normalization in the MC distribution. The best fit values are listed in tab. and the corresponding 
MC fit for the global tracking is shown in fig. ^7|. 

The lateral resolution oi at is derived from the muon resolution parameters a and ji by calcu- 
lating the root mean square of the corresponding distributions of lateral deviations relative to 0, i.e. 
the average lateral distance to the real track. The resulting values for Oi at range from 30 cm for 
global tracking to ^40 cm for OD and ID. This is comparable to the results derived from OPERA 
tracking. The fit value for the neutron mean free path A is close to earlier LVD results I Fill. 

Table 7. Best fit parameters from the numerical fit to the muon-neutron distance (x) distribution. 





o [cm] fx [cm] 


Oiat [cm] 


X [cm] 


* 2 /r.df 


OD tracking 
ID tracking 
GL tracking 


34.4±0.8 20.9±1.0 
34.2±0.8 21.0±0.9 
18.7±1.0 23.0±0.6 


40.3±1.2 
40.2±1.2 
30.0±1.0 


77.4±0.9 
73.9±0.9 
76.6±0.8 


21.4/16 
25.9/16 
31.9/16 
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5. Muon-Induced Neutrons 



Cosmic muons crossing the SSS can produce neutrons through various spallation processes IfRll 
on carbon nuclei. The scattering and thermalization of the neutron is invisible in Borexino as it 
cannot be disentangled from the much larger light emission of the muon. However, the neutrons 
are subsequently captured in the detector with the emission of a 2.2 MeV y-ray if captured on 
hydrogen H or a 4.9 MeV y-ray if captured on carbon C. The typical capture time is of the order 
of 250 (as. Detection of such y-rays with the highest possible efficiency is of primary importance 
for Borexino both for antineutrino detection ([U ED) and for the rejection of cosmogenic isotopes 
produced in association with the neutrons themselves. 

5.1 Neutron detection principles 




time ftis] 

Figure 28. Hit time distribution for a sample IVji event (blue) and the associated neutron detection gate 
(red). 53 neutron captures (spikes) were identified. The smaller graph is a close-up view of the first 30 {is. 

The difficulty for neutron detection arises from the muon after-pulse hits that are generated in 
the PMTs for several tens of |lis after the muon has already left the detector. This is particularly true 
for IVji's, while it is a secondary problem for Bji's. The number of neutrons generated and captured 
in the active volume of the detector after a single muon event can extend up to ^200, introducing 
further difficulty. The challenging task of recording the highest possible number of y-rays within 
a few tens of jus is met as follows. Defining muon event an ordinary ID event flagged also by the 
OD (MTF flag), we programmed the trigger to generate after each muon event a special neutron 
event whose gate is 1.6 ms, i.e. 100 times longer than an ordinary event. Although the muon and 
neutron events are meant to be consecutive a tiny gap could not be avoided and has been measured 
to be (150 ± 50) ns. Hits recorded in the neutron event include the tail of the muon after-pulse hits 
distribution as well as the hits due to neutron captures if present. A sample pair of events where 
a high number of neutron captures is present can be seen in fig. ^3|. A dedicated algorithm has 
been developed to search and isolate cluster of hits due to y-rays that stand out above the "noise" 
background of after-pulse hits. 

Physical events not correlated with the muon may accidentally fall inside the long neutron 
gate and be identified as clusters by the algorythm. They are generally due to 14 C decays corre- 
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sponding to about 50-100 hits and can be confused with partially contained y-rays. However, such 
clusters generally occur in single-cluster events and are mostly rejected by a proper cut introducing 
a negligible loss of signal. 

To test the neutron detection efficiency, a dedicated hardware system has been developed and 
installed: a 500 MHz waveform digitizer (CAEN v 173 1) is fed with the analog sum of all ID PMTs 
properly attenuated. A test performed based on a few months of data returns a neutron detection 
efficiency of ^99% for neutrons produced by IDji's. 



5.2 Neutron capture spectrum 

The energy spectrum for neutron capture y-rays inside the fiducial volume FV (the innermost vol- 
ume, a sphere of radius of 3.1m) is shown in fig. The dashed line is the distribution of the 
number of hit PMTs, n^, normalized to 2000 live channels and shows two prominent features 
around 750 and 1250, due to y-rays from capture on H and C respectively. The solid line is the 
same quantity corrected taking into account the binomial probability of multiple hits: 

AUr = -2000 -In (\ - " h ) (5.1) 

V Jne^vPMTs / 

where M>mts is the number of working PMTs in the run and / ne is the fraction of non saturated 
electronic boards in the specific event. This last correction is necessary because muons that cross 
several meters of scintillator can generate enough light to flood the memory of a fraction of the 
digital electronic boards, which then remain unavailable during the following neutron detection 
gate. However, the number of such boards is recorded and the effect can be corrected as in eq. |5TT| . 

Moreover, neutrons whose production and capture point lie on different sides of the inner 
vessel (either way) can result in partially contained y-rays. For this reason the distribution of 
detected hits for neutron capture y-rays extends down to below the detection threshold visible in 
fig. The structure around 200 hit PMTs is mostly due to neutrons occurring in events with 
many neutrons, induced by showering muons; the neutron visible energy is suppressed by the 
mechanisms described in sec. 1531 . 
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Figure 29. Neutron capture spectrum (see sec. pT2[). Figure 30. Neutron radial distribution (see sec. |53l). 
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5.3 Neutron position reconstruction 

The position reconstruction algorithm works well for neutron clusters belonging to muon events, 
that are not showering and generate normally a lower number of neutrons in the detector. This 
is the case for both Bjli's and IVji's. In general only a slight systematic shift of the reconstructed 
neutron positions to the center of the detector is observed. 

On the contrary, muon events, that are showering and that generate a large number of neutrons, 
produce also a large amount of scintillation light. This leads to two effects, which deteriorate the 
efficiency of the position reconstruction algorithm. 

First, large afterpulses are generated. In the time profile of the signal (fig. ^S|), the neutron cap- 
ture peaks are superimposed on the tail of the afterpulse hits, affecting the position reconstruction 
algorithm as there is no way to discriminate y-ray hits from after-pulse hits. 

The second effect is that a fraction of the ID electronics boards can be saturated and become 
unavailable for neutron detection (as mentioned before). Since PMTs were cabled in bundles, the 
saturated boards correspond to PMTs preferentially located on the half detector closer to the muon 
track. This decrease of detector uniformity aggravates the effect of afterpulse hits on the position 
reconstruction. To overcome this problem and to improve the position reconstruction of neutrons, 
we have developed a new custom firmware for the digital electronics that accepts only 25 hits per 
channel in the muon event, thus preventing the saturation of the memory in the neutron event. Data 
acquired from middle 2010 onwards will benefit from this improvement. 

Fig. shows the radial distribution of neutron capture events. Here we have selected only 
neutron captures past the first 300 |lis after the muon, to limit the influence of after-pulse, and above 
an energy threshold of 200 hits (or lower depending on the level of empty boards), to exclude 
accidental coincidences due to 14 C. The distribution was fitted by the convolution of a volumet- 
ric function with a Rayleigh function that describes the detector resolution. The best fit returns 
# 2 /ndf = 364/191. The goodness of the fit is limited by the fact that a higher energy threshold 
progressively suppresses neutron captures occurring close to the IV surface whose y-rays are only 
partially contained and can introduce distortions in the radial profile at large radii. 

5.4 Neutron capture time 




time [jjls] time [|xs] 

Figure 31. Neutron capture time from cosmogenic neutrons (left panel) and from the Am-Be source 
campaign (right panel). 
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Neutrons produced in muon interactions are slowed down to sub-eV energies in a few scatter- 
ing interactions and this typically occurs within a few tens of ns. Thereafter, the thermal neutron 
can wander in the detector for a relatively long time before capture on H or C occurs. The cap- 
ture time profile shows an exponential behavior (fig. [H], left panel) and depends mainly on the H 
density of the medium. The fit is performed based on an exponential plus a constant term in the 
[100,1500] |is range where the neutron detection efficiency is believed to be uniform. The fit is ex- 
cellent (# 2 /ndf = 134/137) and returns T = (257.4 ±2.4) \is and a background component in the 
sample of about 0.8% (mostly 14 C). However the results are not very stable against radial selection 
of the active volume and against an energy cut, due to the combined effects of the muon after pulse. 
A conservative estimation of the systematical error is about ±8 jas. 

A far more accurate measurement of the neutron capture time has been performed by studying 
an Am-Be neutron calibration source which was inserted into the detector and deployed in more 
than 20 positions within the Inner Vessel [E3I1 . The capture time was determined by the selection 
of coincidences between prompt and delayed signals from the source. The maximal distance of 
correlated events was set to 2 m in order to contain most of the emitted neutrons. The fit of the time 
profile of the source in the center of the detector (fig. [31], right panel) returns T = (254.3 ± 0.8) (as 
(# 2 /ndf = 138/137). The mean neutron capture time obtained from several source positions with 
up to 4.0 m distance from center of the detector is T = (254.5 ± 1.8) jus. 

6. Conclusion 

A key element for the underground Borexino experiment is the suppression of the background 
induced by the residual cosmic muons. A high-efficiency veto and geometrical reconstruction of 
these muons is mandatory for an identification of cosmogenic backgrounds based on time and space 
coincidences with the parent muon. 

The Borexino detector is formed by an inner detector which contains liquid scintillator, the 
central target volume surrounded by a buffer. It is also equipped with an outer detector, a water 
tank, which serves both as an additional passive shielding against external radiation and as an 
active veto, identifying crossing muons via their Cerenkov light emission. The inner detector has 
been largely described in [@j]. The components and characteristics of the outer detector have been 
detailed in the present paper. 

The identification of muons is based on the complementarity between outer and inner de- 
tectors. Different methods of muon identification have been studied and their veto efficiencies 
calculated. The overall Borexino muon identification efficiency reaches 99.992 %. 

The reconstruction algorithms for muon tracks have been described and their performances 
studied using different methods. The resolutions achieved are at the level of ^3.5° — 5° angular 
and ^35-50 cm lateral. 

The cosmogenic neutrons generated by muons have been studied and special techniques have 
been developed to increase the efficiency for vetoing both neutrons and subsequent decays of cos- 
mogenic radioisotopes. 

The study of cosmogenic isotopes induced by muons will be performed in a forthcoming paper. 
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A. The CNGS neutrino beam in Borexino 

The CERN to Gran Sasso (CNGS) neutrino beam [0, EID produces muons in the rock upstream 
of Borexino via charged current interactions of the v^s. Further charged and neutral current 
interactions occur inside the detector. These events are a background to neutrino analysis and 
must be removed from the data set. At the same time they constitute a sample of muons of known 
direction and we use them to determine Borexino muon tagging efficiency as explained in sec. |3.2.1| . 
In addition, since Borexino is just in front of the OPERA detector in Hall C, we exploited the on- 
time events common to the two experiments to test the Borexino tracking algorithms for horizontal 
tracks (sec. fl~2TD . 

The mean neutrino energy is about 18 GeV and the mean energy of the induced muons is about 
10 GeV 1 12511. The expected number of muons per charged current neutrino interaction is one, with 
only 0.3 % of the interactions producing more than one muon. At the LNGS, the neutrino beam 
arrives with an incline of 3.2° below the horizon, corresponding to a a zenith angle of 93.2°; the 
azimuth angle with respect to Hall B axis (parallel to Hall C, where Borexino is located) has been 
measured to be 0.18° fEHl : its width is more than one kilometer. 

For synchronization purposes, both CERN and LNGS have timing system based on GPS, with 
comparable resolution (< 100 ns) and accurately intercalibrated lETIl. The time delays between the 
LNGS master GPS receiver located above ground and the slave units located at each underground 
experiment have been accurately measured. In Borexino the correlation with the beam extraction 
windows can be obtained by comparing the GPS time of each trigger with the ones obtained from 
the CERN database. The time-of-flight distribution shown in Fig. ^ (left) features a clear peak at 
2.4 ms, the expected value for the distance from CERN to LNGS. Fig. (right) shows the same 
quantity on a finer scale allowing to see the ^10 (as spread that corresponds to the extraction time 
duration. The events in the time window [2.40, 2.416] ms are flagged as "CNGS on-time events". 

The fraction of accidental coincidences falling in this window - mostly due to 14 C decays in 
the scintillator - can be easily calculated from the average trigger rate, the number of on-time events 
and the number of extractions during Borexino uptime derived from known beam parameters. We 
obtain / ~1.9 % (3.6 %) if all (only ID) events are considered. This reduces to below ^0.06 % if 
only events above 80 hits are considered as in sec. |3.2.1| . 
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Figure 32. CERN-Borexino time-of-flight from clock time differences. 

Tab. |8| summarizes for each year of CNGS operation the beam intensity and the number of 
on-time events in Borexino. Since some of our calibration campaigns took place during beam time, 
the integrated intensity in the Borexino data sample is generally lower than the intensity delivered 
from CERN. 

Table 8. CNGS beam operation periods; Integrated beam intensity delivered by CERN; integrated beam 
intensity in the Borexino data sample; on-time events detected by Borexino; on-time ID events. Intensities 
are given in units of 10 19 p.o.t. (protons on target). 



year 


beam period 


/ delivered 


/ in BX data 


on-time 


ID on-time 








10 19 p.o.t. 


10 19 p.o.t. 


events 


events 


2007 


Sep 22 


- Oct 20 


0.08 


0.07 


873 


438 


2008 


Jun 18 


- Nov 03 


1.78 


1.33 


17687 


9362 


2009 


May 27 


- Nov 23 


3.52 


2.30 


33206 


17765 


2010 


Apr 20 


- Nov 22 


4.04 


3.16 


45774 


24749 


total 


9.42 


6.86 


97540 


52314 
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